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Abstract 

A (microscopic) static elastoplastic field theory of dislocations with 
moment and force stresses is considered. The relationship between the 
moment stress and the Nye tensor is used for the dislocation Lagrangian. 
We discuss the stress field of an infinitely long screw dislocation in a 
cylinder, a dipole of screw dislocations and a coaxial screw dislocation in 
a finite cylinder. The stress fields have no singularities in the dislocation 
core and they are modified in the core due to the presence of localized 
moment stress. Additionally, we calculated the elastoplastic energies for 
the screw dislocation in a cylinder and the coaxial screw dislocation. For 
the coaxial screw dislocation we find a modified formula for the so-called 
Eshelby twist which depends on a specific intrinsic material length. 

Keywords: elastoplasticity, dislocations, moment stress 
PACS: 61.72.Lk, 62.20.-x, 81.40.Jj 



1 Introduction 

The traditional description of the elastic field produced by a dislocation is 
based on the classical theory of linear elasticity [1-4]. This approach works 
for the strain and stress field far from the core quite well. However, the 
components of these fields are singular at the dislocation line and this theory, 
often applied to practical problems, misses the important feature of plasticity. 
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This is unfortunate since the dislocation core is the most important region in 
metallic plasticity. 

Therefore, it is quite natural to think of dislocation theory as a field theory 

of elastoplasticity [5-11]. In this framework, it is possible to understand the 
dislocation density as the elastoplastic field strength like the magnetic field 
strength in Maxwell's theory. The presence of dislocations leads to a specific 
response with the dimension of a moment stress, so that the excitation or 
response quantity in elastoplasticity can be identified with a localized moment 
stress tensor. In order to obtain a physical field theory of dislocations, we have 
to choose a constitutive law between dislocation density and moment stress 
which is compatible with the symmetric force stress. We use the constitutive 
law for the moment stress which is proportional to the Nye tensor. A new 
material constant enters the constitutive law. Additionally, it is possible to 
define a characteristic internal material length in this theory by means of the 
new material constant. 

On the other hand, there are other non-standard continuum models of 
dislocations as, e.g., the nonlocal continuum model [12-15] and strain gradient 
elasticity [16-21]. In these approaches also a characteristic length (gradient 
coefficient or nonlocality parameter) appears. In the case of a straight screw 
dislocation in an infinite medium, the force stresses in nonlocal elasticity [12- 
15] and in gradient elasticity [19, 20] agree with the gauge theoretical ones [10, 
11, 22, 23]. These solutions have no singularities at the dislocation line. 

The plan of the paper is as follows. In Section 2, the linear elastoplas- 
tic field theory of dislocations [9, 10] is recalled. In this paper we use the 
widespread formalism of tensor analysis because we believe that this formula- 
tion of the elastoplastic field theory is comprehensible and useful for material 
scientists, experimental physicists and engineers. We give the relation of this 
elastoplastic theory to Kroner's dislocation theory [24,25]. In Section 3, we 
investigate screw dislocations with different boundary conditions in this frame- 
work. 

2 Linear elastoplastic field theory of dislocations 

In the elastoplastic field theory, the incompatible distortion tensor 

Bij = dj^i + (pij, (1) 

plays a fundamental role [9, 10]. Here is the mapping function from the ideal 
or reference configuration to the final or deformed configuration. The first part 
of Eq. (1) describes the integrable part of the distortion, while (f)ij is, in general, 
the nonintegrable distortion which can be identified with the gauge potential of 
dislocations in the framework of the r(3)-gauge theory. The form of Eq. (1) 
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can be understood as a minimal replacement of the compatible distortion, 
Bij = dj^i, in T(3)-gauge theory and is obtained from the translational part 
of the generalized affine connection [9, 10]. The distortion Bij describes the 
local distortion from the ideal or reference state to the dislocated state which 
is now anholonomic. 

For simplicity, we will consider here a linear isotropic theory of dislocations. 
In this case, we do not distinguish between covariant and contravariant indices. 
A nonlinear clastoplastic field theory of dislocations has been proposed in 
Rcfs. [9, 10] (sec also [11]). 

If we use the weak field approximation (linearization) 

d — Xi -\- Ui, Bij = dij + Pij, (2) 

where Ui is the displacement field, the linear elastic strain tensor is given by 
means of the incompatible distortion tensor as 

Eij = P^ij) = ^ {diUj + djUi + (f)ij + ^ji) . (3) 
The elastic strain energy of an isotropic material reads 



W = - /J, (^SikSji + SiiSjk + -^_2j^ ^ij^ki^ EijEki, 



(4) 



where ^ is the shear modulus and u Poisson's ratio. The (symmetric) force 
stress is the response quantity (excitation) with respect to the strain and is 
given by 

The total strain Ej^ and distortion Bfi are defined in terms of the total dis- 
placement field uj. In the presence of dislocations, the total strain is not 
completely elastic, but a part of it is stress-free or plastic [26] so that 

Efj ^ d^iuj^ = E,, + Efj, (3fj = djuf = (3i, + Pfj. (6) 

The skew-symmetric part of the distortion tensor defines the elastic rotation 
of a dislocation [26] 

= "2 ^ijkPjk- (7) 

The elastoplastic field strength is defined by means of the distortion tensor 
(gauge potential) as follows 



Tijk '■ — djPik dkf^ij — 9j(f)ii~ dk^'ij, T^jk — Tikji 



(8) 
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and 

Tijk = -djf3g + dkf3tj, (9) 

where Pfj is the plastic distortion tensor which may be identified with —(pij. 
The usual dislocation density aij, which was originally introduced by Nye [27], 
is recovered by 

"ij •= 2 ^ ^jkidkPii = -ejkidkPu- (10) 

In this approach the dislocation density may not assumed as a delta function. 
It satisfies the translational Bianchi identity 

djaij = 0, (11) 

which means that dislocations cannot end inside the body. Note that the aij 
defined here is the transpose of the dislocation density aji used by Kroner [24] 
(and —aji in [25]). 

We use the Lagrangian of dislocations in the following form [9, 10] 

■^disl ~ '^Tiji; Hiji;, (12) 

with the moment stress i7y ^ as the elastoplastic excitation (response quantity) 
to the dislocation density (torsion)^ 

Hijk = -2^7j; — = — {Tijk - Tjki - Tkij - 2SijTiik - 25ikTiji) . (13) 
olijk 2, 

Here, Hijk = —Hikj- The coefficient ai has the dimension of a force. If we use 
the relations iJj^ = \ejkmHijk and Tijk = ^jkmO-im-, we find from Eq. (13) the 
relation between the moment stress tensor Hij and the Nye [27] tensor Kij as 

Hij = ai (^ji - ^ dijakk^ = aiKij- (14) 

Hence, this moment stress tensor Hij is proportional to the Nye tensor (see 
also [28]). The diagonal components Ha describe twisting and the nondiagonal 
components of Hij bending moments localized in the dislocation core region. 
Alternatively, the dislocation Lagrangian (12) can be written as 

Aisi = - Y aijKij. (15) 



^Here, the so-called Einstein-choice is used (see [10]). 
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In the framework of elastoplastic field theory, we require that the modified 
stress field of a dislocation has the following properties: (i) the stress field 
should have no singularity at r = 0, (ii) the far field stress ought to be the 
stress field of a Volterra dislocation (background stress) aij which satisfies the 
condition djaij = 0. The boundary Lagrangian is given by a so-called null 
Lagrangian [29] and by using the minimal replacement djUi Pij: 

W\yg = dj{aijUi) = {djcrij)ui + a^jdjUi -> aijjSij. (16) 

The Euler-Lagrange equations for £, = Cdisi — W + VFbg take the form 

— (force equilibrium), (17) 

dkHikj = o'ij (moment equilibrium), (18) 



where the so-called effective stress tensor. 



(19) 



is the driving "source" for dislocations. The left hand side of Eq. (18) reads 
explicitly 

= AEij - {djdkEik + didkEkj) + SijdkdiEki + didjEkk - SijAEkk 
= (inc-E)-^- = -eikiejmndkdmEin, (20) 

and 



dk%,]k = 0. (21) 

Consequently, the antisymmetric moment equilibrium or equilibrium of couple 
stress (21) is trivially satisfied in the Einstein choice and, thus, no antisymmet- 
ric force stresses appear. Additionally, we obtain the equation for the elastic 
strain Eij as 

ai{AEij — {djdkEik + didkEkj) + SijdkdiEki + didjEkk — ^ijAEkk} = ^ij- 

(22) 

This is the incompatibility condition of the elastoplastic theory of dislocations. 
By means of the inverse of Hooke's law 
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and the equilibrium condition djaij = 0, we obtain the field equation for the 
stress field 

AcTjj + {didj - SijA)akk = ^ij, = (24) 

For aij = 0, this equation is the Beltrami equation. The coefficient k has the 
dimension of a reciprocal length. Now, we compare Eq. (22) with Kroner's 
incompatibility condition [24, 25] 

— ^ikl^jmndkdmEln = Vij^ (25) 

where the so-called incompatibility tensor r]ij is defined by 

Vij ■= -^{^ikldkoiij + ejkidkOiii) ■ (26) 

Then we can identify 

riij = —aij. (27) 
ai 

Therefore, Eq. (22) represents the proper gauge theoretical formulation of 
Kroner's incompatibility equation. However, Edelen's field equation (see, e.g., 
Eq. (9-8.7) in Ref. [29]) 

Apij-djdil3u = —^ij, (28) 

0-1 

cannot be interpreted in this way, in contrast to his claim. 
We may define the tensor 

Qij '■= —^jkidki^ii (29) 

which can be identified with the linearized Einstein tensor in a Riemannian 
space. Here, Qij is determined by the dislocation density. Using Hijk = ejkiHu 
and Eq. (14), the field equation (18) has the shape of a linearized Einstein-type 
field equation according to (see also [9-11]) 

= ^ij with djaij = 0, dj&ij = 0. (30) 

3 Screw dislocations in the elastoplastic field theory 

3.1 Screw dislocation in a cylinder 

We now consider an infinitely long screw dislocation whose dislocation line 
and Burgers vector coincide with the axis of an isotropic cylinder of radius R. 
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Due to the symmetry of the problem, we choose the Burgers vector and the 
dislocation line in z-direction: bx = by = 0, bz = b. Using the stress function 
ansatz [4, 24, 25] , the elastic stress of a Volterra screw dislocation, in Cartesian 
coordinates, reads 

(^xz = CFzx = -Oy^ = (^yz = <^zy = Ox^ ^ 27r ^' ^ ' 

or, in cylindrical coordinates, 

O O r\ X. /n<->\ 

(^zip = (^<pz = Or^ = (32) 

where = x'^ + y'^. Here, $ is the well-known stress function of elastic torsion, 
sometimes called Prandtl's stress function. It is given by (see, e.g., [25]) 

$ = ^lnr. (33) 

This is Green's function of the 2-dimensional potential equation 

A$ = /x6(5(r). (34) 

Obviously, the "classical" stress fields are singular at the dislocation line. For 
simplicity, we assume that the cylindrical surface r = i? is free of tractions, 
i.e.. 



o o 



= for r = i?. (35) 
For the modified stresses we make the ansatz 

-dyF, ayz = CFzy = d^F, (36) 



'xz — zx 



where F is called the modified Prandtl stress function. We turn to Eq. (24) 
and use akk = 0. In this way, the field equation for the force stress of a linear 
screw dislocation reduces to an inhomogeneous Helmholtz equation, 

(l - K~'^A) Gij = aij, (37) 

with the "source" term given in terms of the well-known solutions of classical 
dislocation theory for the same traction boundary-value problem. Thus, for 
the elastic strain fields, 

(1 - K-^A) Eij = Eij, (38) 
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where Eij is the "classical" strain tensor. Substituting the stress functions 
into (37), we get 

(l-K-2A)F = ^lnr. (39) 

The solution of the modified stress function of a screw dislocation is given by 

F = ^{lnr + KoiKr)}, (40) 

where Kn is the modified Bessel function of the second kind and of order n. 
Consequently, we find the force stresses 

(^xz = cTzx = ^{^~ KrKi{Kr)Y Gyz = ^^y = ^ ~ nrKi^Kr)^ 

(41) 

and in cylindrical coordinates 

crz<p = cr<pz = i^rKi (kv) | . (42) 

If we choose the coefficient k as (a = lattice constant) 

^ 0.4a, (43) 

the stress azy has its maximum 0.4K/Lt6/27r = /x6/27ra at x ~ 1.1/ k = 0.44a 
whereas the minimum — 0.4K/^6/27r = —nb/2TTa lies at x ~ —1.1/k = —0.44a. 
This is in agreement with the Peierls-Nabarro model. Note that the maximum 
stress may be identified with the theoretical shear strength. The core radius 
Tc is then given as Tc c± 2.4a (Figure 1). This estimate of k is in conformity 
with the result in Eringen's nonlocal elasticity theory. He pointed out that for 
the choice of ~ 0.4a one has an excellent fit between atomic dispersion 
curve and the nonlocal result [12-15]. However, the factor k should be fitted 
by comparing predictions of the theory with experimental results. In general, 
the parameter k can be used to determine the width of a dislocation and the 
amplitude of the force stress. 

Let us now calculate the distortion of a screw dislocation. The distortion 
Pij is given in terms of the stress function 

—dyF + 2iJLUJi 

Ai = £ I - 2a«^2 I , (44) 
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Figure 1: Modified stress field ayz{x,0)/{iJ,b/2Tr) (solid), with 2± 0.4, 
a = 1, and classical stress (dashed) of a screw dislocation plotted over the 
x-coordinate. 

where uji and L02 are used to express the antisymmetric part of the distortion. 
Eventually, coi and 0^2 are determined from the conditions 

oiyx = Tyyz = -^dy (2/xa;2 - d^F) = 0. (46) 
One finds for the distortion tensor of the screw dislocation 

and for elastic rotation vector 

h X ( "1 ^ y ( 1 

a;^ = W2 = ^ - KrKi{Kr)j, Uy = = — - KrKi{Kr)y 

(48) 

The rotation vector is in agreement with the result calculated in the Cosserat 
theory [30]. 

By means of the distortion tensor it is possible to determine the effective 
Burgers vector as 

bz (r) = / [Pzxdx + P,ydy) = b{l - nrKi {kv)}, (49) 

where 7 is the Burgers circuit. This effective Burgers vector differs appreciably 
from the constant value b in the core region from r = up to r ~ 6/k. Thus, 
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it is suggestive to take Tc ^ 6/k as the characteristic length (dislocation core 
radius). Outside this core region, the Burgers vector reaches its constant value. 
Accordingly, the classical and the elastoplastic solution coincide outside the 
core region. 

If we use the decomposition (6) of the elastic distortion into the total and 
the plastic distortion, the total displacement field turns out to be 



= (fi^l — KrKi{Kr)^ , = arctan—. (50) 



Thus, the proper incompatible part of the total distortion is the plastic dis- 
tortion which is confined in the dislocation core region 

l^zx = -^ x<pKo{Kr), (5^y = — yipKo{Kr). (51) 

It fulfills Eq. (10). When 0, the plastic distortion (51) converts to 

= bxipS{r) and = byipS{r). In our framework, Eq. (50) describes the 
"atomic" arrangement in the dislocation core region in contrast to classical 
elasticity which leads to an unphysical singularity. The incompatible elastic 
and plastic strains read, respectively. 



■xz - ^zx = ^|l - KrKi(Kr)Y Eyz = E^y = ^ ^|l - KrKi{Kr)^, 



Etz — E. 

47r L - - j - -~ ~« L 

(52) 



and 



5^2 u 2 

eL = C = 4^ x^Koinr), = = — y^Mi^r). (53) 



Now we are able to calculate the dislocation density by means of the dis- 
tortion tensor. We obtain 

1 

azz = Tzxy = -T,y^ = -AF = — Koinr). (54) 

/J, ZTT 

Of course, this dislocation density satisfies the conditions (11), (26) and (27). 
In the limit as 0, the elastoplastical result (54) converts to the classical 

dislocation density azz = b6{r). The localised moment stresses caused by the 
screw dislocation can be expressed in terms of the dislocation density as 

FT- - - ii F7 - - _ii H- - 

J^zz — (^If^zz — 2 ^^^^ ^xx — G-lf^xx — 2 ^^^^ ^yy — ^l^yy — 2 ^^Z"! 

tx/ tv tv 

(55) 

and 

Hkk = --^ oizz, Kkk = -7; oizz- (56) 
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Accordingly, moment stresses of twisting-type occur at all positions where the 
dislocation density azz is non-vanishing. When — > 0, the moment stresses 
vanish. 

We next define the characteristic intrinsic material length by 

Rc measures the region over which the dislocation density (torsion), the Nye 
tensor and the moment stress are appreciably different from zero. It can be 
called the plastic penetration depth. Then i?c is the internal length which 
measures the plastic region where (ifj ^ 0. In this region the local atomic 
configuration is fundamentally different from that of the defect-free part of 
the crystal. An important consequence is that the mechanical but also the 
electrical properties arc different from those in the undefected crystal. 
The strain energy density produced by the screw dislocation is 

W = \ aijEij = ^ (1 - KrK,{Kr))\ (58) 

We find the strain energy stored per unit length of the screw dislocation in 
the cylinder as 

rR j-2-K 
EstraXn = drr difiW 

Jo JO 

47r 

where R is the outer cut-off radius. We use the limiting expressions for r — > 0, 



2 2 

■{inr + 2Ko(Kr) + ^-{Ki{Krf - Ko{Kr)K2{Kr))j^^ , (59) 



i^o(/^r)«.-7-lny, K,{nr)^^, K2(«r) « + (60) 

where 7 = 0.57721566 ... is the Euler constant. The final result reads 
fib"^ (, kR 1 



^strain = ^{^^^+^-2 + 

[KiinRf - Ko{KR)K2{KR))y (61) 



K^R^ 

+ — — 



Thus, we obtain a strain energy expression which is not singular at the dislo- 
cation line and at R. 

The dislocation core energy density produced by the screw dislocation is 
given by 

Vcore = y OLzzi^zz = -^o(Kr) . (62) 
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The dislocation core energy per unit length of the screw dislocation is 

rR r2iT 
-Score = drr dip Vcove 

Jo Jo 



+ K^R^{Ko{KRf - Ki{KRf)y (63) 



Stt 

_ 

~ Svr 

Eventually, the total energy per unit length of the screw dislocation in a cylin- 
der is given by 

-^screw ~ -E^strain "I" -E'core 

= ^ |ln^ + 7 + 2Ko{kR) - kRKq{kR)Ki{kR) \ . (64) 
47r [2 J 

It depends on the radius R of the cylinder. 

3.2 A dipole and parallel screw dislocations 

Because the strain and stress fields in our approach have not any singular- 
ities, it is interesting to discuss the stress fields of two (anti)-parallel screw 
dislocations {A) and [B) near the dislocation cores. These two dislocations 
lie in the (xy)-plane at x = 0, y = {A) and x = d, y = (B) with the 
Burgers vector extending in z-direction. Eventually, the stress fields of these 
two screw dislocations are given by the superpositions of the corresponding 
fields of individual ones as 

o-L = {~^{^ ~ KrAKi^KrA}^ + ^^{l - KrB-f^i(KrB)|^ . (65) 

Here r\ = + and = {x — df' + y^. For 6^ = hs and 6^ = —6b, this 
configuration describes two parallel screw dislocations and a dipole of screw 
dislocations, respectively. Again, in contrast to the classical stress fields, the 
modified stress fields dislpay a smooth shape and have no singularity (Fig- 
ure 2). 

We may discuss the stress field of a wall of screw dislocations. A finite wall 
of (anti) parallel screw dislocations lying in the plane x = is obtained by a 
superposition of the individual screw dislocations. The stress field of this wall 
is given by 



a. 



sr^ ^J'On X ( \ II s:^ pLbnyn 



(66) 



n=0 " n=0 " 
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(b) 

Figure 2: The stress component alz{x,0) /{fibA/^Tr), with d = 5, = 0.4 
and a = 1 (solid), (a) of a pair of parallel screw dislocations, (b) of a screw 
dislocation dipole. The dashed curves represent the classical stress component. 

where AT + 1 is the number of the dislocations in the wall, bn is the Burgers 
vector of the individual dislocations, = + y^, yn = y — nd and d is the 
spacing of the dislocations. 

Because the stress fields have no singularities, it is possible to study the 
behaviour and the interaction between the screw dislocations in this frame- 
work. 

3.3 A coeixial screw dislocation in a cylinder 

Now we consider a screw dislocation in a cylinder of finite length in z-direction. 
Up to now, this problem was only discussed in classical elasticity (see [2, 4, 31- 
33]). We use the solution (42) which fulfills the boundary condition at the 
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cylindrical surface r = R. On the end of the surfaces the solution (42) produces 
the torque 

Mz = dip / drr a^^ 
Jo Jo 

HbR^ r 4 



{i-^^r^ + ^^^M)} (67) 



on the upper basis and —M^ on the lower basis of the cylinder. Thus, Eq. (42) 
is not the correct solution for a finite rod and it has to be modified. The torque 
is removed by adding an equal and opposite one by means of an additional 
background condition 

Ux = '^yZ, Uy = -TXZ, 

^'xz = ^'zx = 'Tl^y, ^'yz = ^'zy = -TfiX, (68) 

with a coefficient r which is the so-called Eshelby twist. Then the correct 
boundary condition for the stress field is 

o o l^b y o o fib X 

(^xz = <^zx = + rny, (Tyz = = ~ '^i^^- y^^j 

Eventually we find from Eq. (37) 

cTxz = CTzx = - KrKi{Kr)^ +Tny, 

IJ,b X ( / \1 

o-yz = o-zy = ^ ^ 1 1 - KvKi {Kr) j - TflX, (70) 

and in cylindrical coordinates 



_ Ijb 1 



|l — KrKi{Kr)^ — Tfir. (71) 



' Z(p — '^(fiZ n 

ZTT r 

We determine r from the condition of vanishing torque, = 0, with 

Mz= I d(p drr Uzip 
Jo Jo 

We find the "modified" Eshelby twist 



l-;^ + 2i^2(Ki?)}. (73) 
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The first term is the "classical" Eshelby twist (see [32]) and the last two 
terms are correction terms in our approach. The correction terms depend on 
the intrinsic material length k~^. The Eshelby twist is an observed quantity 
in a thin metal whisker which has a screw dislocation along its axis [4, 34, 
35]. It is important in transmission electron microscopy since it reveals screw 
dislocations normal to the foil under consideration by a characteristic contrast. 
In principle, if the "correction" terms give a measurable contribution, the 
constitutive constant k could be determined in such experiment. 

Finally, the total displacements and stresses for the axial screw dislocation 
in a free rod are given by 



< = -^{'-;^ + ''^^^''^)}' ^^=2^^{^-^rK^{nr)}, (74) 
az^ = a^z = ^l{l- nrK.iKT) - ^ (l - + 2K2{kR))}. (75) 

The rotation vector of the screw dislocation in a whisker reads 

Uz = —Tz. (76) 
If we use the definition of the rotation gradient (de Wit's bend-twist tensor) [26] 

kij = djUJi, (77) 

the far-reaching rotation gradients of a coaxial screw dislocation in the whisker 
read 

kyy = ^{{x^ - y^){l - KrKiiKr)) + KWr^KoiKr)^ + ^T, 
kxy = kyx = -^;^xy|2(l - KrKi{Kr)) - K^r'^KQ{K,r)^, 



kzz — ''') 
1 
2 



kjj = -^Oizz- (78) 



For r = 0, these rotation gradients are in agreement with the expressions 
calculated within the theory of Cosserat media (see [36-38]). Additionally, 
we find the relation Kjj + kjj = 0. Thus, the sum of the Nye tensor and the 
rotation gradient, which may be considered as the total "lattice" curvature or 
total bend-twist tensor, is a deviator, that is, it is traceless. 
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It is interesting to note that the strain energy of an axial screw dislocation 
is influenced from the free-surface terms. Eventually, the strain energy of the 
screw dislocation in a finite cylinder of outer radius R is given by 

^strain =^{^'^^+^-2^ 2Kq{kR) 

K?R^ ^ 
+ -^{KiiKRf - Ko{KR)K2{KR))j 

+ 4ifo(/tii)iC2(Kfl) + --Ki{^R)K2(kR)}. (79) 

The terms in the third and fourth line of Eq. (79) arise from the "modified" 
Eshelby twist. Note that the first term in the third line of Eq. (79) is the 
correction term resulting solely from the "classical" Eshelby twist, see [4]. 
The core energy (63) is not influenced by the Eshelby twist. 

4 Conclusions 

In this paper wc have considered a linear elastoplastic field theory of dislo- 
cations. The total Lagrangian of elastoplasticity has the following symbolic 
form 

jC = C f3'^ + Aoi^ + boundary terms, (80) 

where A, C are material tensors. It is given in terms of the both physical 
state quantities, namely the elastic distortion (3 and the dislocation density 
a. Comparing the dimensions, we note that (3 is dimensionless and a has 
the dimension of an inverse length. Therefore, the second term in Eq. (80) 
contains a factor = l/(length)^ when compared with the first one. In this 
theory, localized moment stress is the specific response quantity to dislocation 
density. Thus, an elastoplastic dislocation theory is a theory of force and 
moment stress. Additionally, we have shown that the field equation of /3 is 
really the gauge theoretical modification of Kroner's incompatibility equation. 

A new material coefficient 04 enters the constitutive relation between dislo- 
cation density and moment stress. It defines a new characteristic length k^^. 
In general, this length specifies the plastic region and has to be estimated 
by means of experimental observations and computer simulations of the core 
region. In the limit ^0, one recovers classical elasticity. 

Exact analytical solutions for the total displacement field, elastic and plas- 
tic strain, stress fields, and total energy of a screw dislocation in a cylinder 
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and of a coaxial screw dislocation have been reported which have not any sin- 
gularity in the dislocation core region. Consequently, it would be interesting 
to compare them with experimental results and related computer simulations. 
For the coaxial screw dislocation, we found a modified formula for the so-called 
Eshelby twist. The results for the screw dislocation in a thin whisker should 
be useful for torsion experiments of thin wires. 

Note that the present theory of elastoplasticity is intimately connected 
with the so-called gradient theory of dislocations [18, 19] and Eringen's nonlo- 
cal elasticity [12-15]. The stress field of a screw dislocation in elastoplasticity 
coincides with the solution in Eringen's nonlocal elasticity [12-15] and in gra- 
dient elasticity [19,20]. The internal characteristic length is thus connected 
to the gradient coefficient and to the nonlocality parameter. Obviously, the 
equations (22) and (37) have the form of gradient theory equations of second 
order. The new material coefficient ai has the meaning similar as a gradient 
coefficient. 
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